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APOLLO MISSION F 
PERFORMANCE ANALYSIS 
OF RENDEZVOUS CHARTS 


1.0 Summary 

A digital analysis of the F Mission backup rendezvous 
charts has been performed to verify their ability to predict 
CSI, CDH, TPI, and midcourse corrections in the presence of 
system and application errors and trajectory dispersions. 
Procedures for data acquisition were as defined in the current 
LM Rendezvous Procedures Document (Reference l). The charts 
simulated were those from the LM 4 onboard data package. The 
mission situation simulated represented the PNGS inoperative, 
rendezvous radar information available from the tapemeter, 
attitude data from the AGS and control by the CES. The study 
showed that the standard deviation of the arrival time at 
TPI was 4.45 minutes. The mean AT from TPI to intercept with 
braking and line-of-sight control was 45.9 minutes with a 
standard deviation of 1.0 minutes. If no braking or line-of- 
sight control was executed after the second midcourse correction, 
the mean miss distance would have been 0.27 NM with a standard 
deviation of .4l NM. The mean total translational AV required 

i 

after insertion was 150.5 ft/sec with a standard deviation of 
21.0 ft/sec. 
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2.0 Discussion 

This analysis was performed to determine the ability 
of the charts to predict all maneuvers after insertion required 
to complete the LM-active F mission rendezvous. One hundred 
independent sets of initial conditions (iC's) were selected by 
adding random errors to the 6 components of the nominal relative 
state vector between the LM and CSM. Forty percent of the 
relative error was incorporated into the inertial state vector 
of the GSM, while the remaining sixty percent of the relative 
error was incorporated into the inertial state vector of the LM. 

The 100 cases were run including system and execution 
errors with braking and line-of-sight (LOS) control. The runs 
were repeated omitting braking and LOS control to obtain miss 
distances. The results of the 100 Monte Carlo runs were pro- 
cessed to obtain statistical data for the parameters of interest. 

The first 50 cases were also run without system or applica- 
tion errors to establish the theoretical chart capabilities. 
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2 » 1 Study Rules 

The following ground rules consistent with F Mission 
procedures and planning were used in this study: 

1. GSI occurs 51 min 01 sec after insertion. 

2. CDH occurs 58 min 01 sec after CSI. 

3. TPI was assumed to occur 9 min after the 

elevation angle of 19.40 degrees was reached. 

4. The first mideourse correction (MCI) occurred 
15 minutes after TPI and the second midcourse 
correction (MC2) occurred at 30 minutes after 
TPI. 

5. CSI and CDH were burned in local vertical 
coordinates using impulsive thrust. 

6. TPI, MCI, MC2, braking, and LOS control were 
executed along and normal to the line-of- sight 
using finite thrust, burning each component 
individually. 

7. 130 degrees of CSM orbit travel between TPI 
and TPF. 

8. No out-of-plane corrections were made prior to 
LOS control during the braking phase. 
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2.2 Di gital Program 

In the analysis, functions describing the backup 
rendezvous charts were programmed into a digital routine 
which integrated the equations of motion of two particles 
about a triaxial lunar gravity model. These functions 
allowed simulation of the procedures for using the backup 
rendezvous charts by incorporating elevation angle (AGS 
address 304), range, and range rate from the tapemetor at 
the times called for by the backup data sequence. System 
errors were included in the data taken at each point and an 
appropriate error of application was included in each maneuver. 

The runs included effects of both bias and random errors 
as defined in Section 2 . 5 . Bias errors were selected by the 
program at the beginning of each run and held constant for that 
run. Random errors were selected at each point where data 


were taken. 
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2,3 Charts 

The charts modeled in the study were those which will he 
carried on the F Mission* and are shown in Figures 1-6. The 
data sequence was obtained from Reference 1 and is summarized 
on the relative motion plot of the nominal trajectory from CSI- 
30 min to TPF (Figure 7). 

2.3.1 CSI 

The CSX chart solution is based on a third order Maclaurin's 
expansion of four variables for the delta V at CSI (range rate 
at 30, 20, and 10 minutes prior to CSI and range 10 minutes prior 
to CSI). The coefficients are determined by the simultaneous 
solution of several expansions, each representing a dispensed tra- 
jectory prior to CSI. 

2.3.2 CDH 

The CDH chart solution is based on a second order Maclaurin*s 
expansion of three variables for the delta V at CDH (range rate at 
36, 23, and 10 minutes prior to CDH). The coefficients are deter- 
mined by the simultaneous solution of several expansions, each 
representing a dispersed trajectory prior to CDH. 

2.3.3 TPI and M CC 

The TPI charts solve for the relative position and velocity 
at TPI resolved into normal and along the line-of-sight (LOS) 
coordinates. The measured relative conditions are differenced 
from the required conditions for intercept in 130 degrees of 
orbit travel. Information required for the TPI charts is the 
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elevation angle at 9 and 5 minutes prior to TPI and range 
and range rate at 5 minutes before TPI. In a similar manner, 
the MCC charts maintain the time of TPF consistent with the 
TPI maneuver. Data are taken for the midcourses at 9 and 13 
minutes after TPI for MCI and at 2k and 28 minutes after TPI 
for MC2. The same measurement sequence as used at TPI is used 
for both midcourse corrections. 
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2.4 Initial Conditions 

The initial conditions were generated by perturbing the 
nominal state vectors of the LM and CSM with errors supplied 
by a relative covariance matrix. Forty percent of the relative 
error was applied to the state vector of the CSM, while the 
remaining sixty percent of the relative error was applied to the 
state vector of the LM. 

The nominal vectors relative to the LM orbital plane were 
derived from Reference 2 and are summarized as follows: 


Table 2-1 
Initial Conditions 
CSI-33 min. 

(103 hrs 0 min 46 sec) 


LM Altitude 

CSM Altitude 

LM Total Velocity 

CSM Total Velocity 

LM Flight Path Angle 

CSM Flight Path Angle 

In Plane Central Angle 

LM Latitude 

Out of Plane Distance 

Out of Plane Velocity 

Heading Angle (relative to equator) 


131859.3 feet 

349583.7 feet 
5458.4 ft/sec 
5350.0 ft/sec 
•935 deg 
.002 deg 
11.869 deg 
0.0 deg 
0.0 feet 
0.0 ft/sec 
0.0 deg 


* 
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The covariance matrix used for initialization of the 
Monte Carlo runs is as follows; 


Table 2-2 
Covariance Matrix 


254189300. 

812877.2 

44433410. 

-53475.3 -424.228 

214069. 

812877.2 

180909.9 

271115. 

-81.066 -2.29299 

644.7976 

44433410. 

271115. 

10541620. 

-6910.21 -113.378 

36791.17 

-53475.3 

-81.066 

-6910.21 

13.45366 .05168073 

-45.6078 

-424.228 

-2.29299 

-113.378 

.05168073 .009133634 

-.343047 

214069. 

644.7969 

36791.17 

-45.6078 -.343047 

I8i.24l4 


It has since been learned that this matrix, obtained from 
Reference 3j reflects trajectory dispersions approximately 4 times 
greater than those which can be expected in the actual mission. 
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2.5 Errors in Sensors and Execution 

The Id errors in sensor and maneuver executions were: 
1. System Errors 


A. 

Noise 



1. Range 

.333$ 


2 . Range Rate 

.333 ft/sec 

B. 

Biases and Drifts (constant 

for a given run) 


1. Initial Pitch Bias 

.1 deg (assumes 
calibrated OOAS) 


2. Pitch Drift Rate 

.18 deg/hr 

Execution Errors 


A. 

Reading Tapemeter 



1. Range Rate 

.25 ft/sec 


2. Range 



a. Outer Scale 

2400 ft 


b. Middle Scale 

IOC ft 


B. Application of Burns .25 ft/sec (per axis) 

C. Time Measurements .5' sec 

D. Elevation Angle .12 deg 

(Error in pointing 

Z-axis at target) 




2-9 


2.6 Braking Schedule and LOS Angular Rate Deadbands 

The braking schedule used in this simulation consisted 
of five gates and a lower limit on the range rate. The first 
gate was at 13500 feet. At this point only LOS control was 
executed because the allowed range rate was 80 ft/sec. The 
second gate was at 6000 ft with an allowed range rate of 30 
ft/sec. The nominal range rate at this range was 32 ft/sec. 
The remaining gates were 20 ft/sec at 3000 ft, 10 ft/sec at 
1500 ft and 5 ft/sec at 500 ft. The lower range rate limit 
consisted of a straight line connecting 20 ft/sec at 13500 ft 
and 0 ft/sec at intercept. 

Both the upper and lower range rate limits are shown in 
Figure 8. LOS control procedures were simulated by sampling 
inertial drift of the LOS inplane and normal to the orbit 
plane every 15 seconds beginning at a range of 13500 ft. When 
the LOS rates exceeded .3 mr/sec at a sampling time, thrust 
was applied in the appropriate axis in increments of 1 second 
until the LOS rate was reduced below the threshold. The 15 
seconds were allowed to elapse before sampling again. 
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3 . 0 R esults 

Several sets of Monte Carlo runs were made to obtain 
statistical data for determination of the effects of errors, 
trajectory dispersions, and braking on the size of maneuvers, 
arrival time at TPI, and total translational AV. The sets 
of runs are identified in the following table: 

Table 3-1 
Run Summary 


Set 

Number 
of Runs 

Errors in 

Maneuver 

Solutions 

Maneuvers 

Applied 

LOS Control 
and Braking 

A 

100 

Yes 

With Errors 

Yes 



No 



B 

100 

Yes 

With Errors 

No 

C 

50 

No 

No Errors 

No 

D 

50 

No 

No Errors 

Yes 


Solutions 

for the maneuvers 

in Set A were 

obtained with 


errors. In addition, conic solutions were obtained for each 
maneuver except CSI. These were used to establish nominal 
AV values, which were used as a basis to compute the A AV 
values in Table 3-3. In the case of CSI, a chart solution with 
no sensor and reading errors was used for this purpose. All 
solutions which were actually applied in Set A, however, were 
arrived at using sensor and reading errors. The runs for Set B 
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ware identical to Set A, tut with braking and LOS control 
omitted to establish miss distance. 

Sets C and D were run to establish baseline data for 
chart performance. It was felt that a reduced number of runs 
would suffice to obtain statistically meaningful results 
since only initial conditions were varied. Examination of 
significant parameters such as maneuver solution and miss 
distances revealed nearly normal statistical distributions, 
confirming the adequacy of 50 runs for those sets. The runs 
for Sets 0 and D were made with the same initial conditions as 
the first 50 runs of Set A. 
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3.1 Maneuver Values 

The nominal solution, along with the average, mean, and 
standard deviation for each maneuver in Set A and Set D are 
shown in Table 3-2 on Page 3-4. The data given for Set A are 
the solutions with errors. The average, mean, and standard 
deviation for the difference between the error solution and 
the theoretical solution computed for each maneuver in Set A 
are shown in Table 3-3 on Page 3-5 • 

It can be noted from the data in Table 3-3 that the chart 
solutions with errors progressively decrease in accuracy from 
CSX to CDH to TPI. The trend then reverses with MCC1 more 
accurate than TPI and MCC2 more accurate than MCC1. 

In addition, it should also be noted that the value listed 
for the nominal C8I solution in Table 3-2 differs from the 
current value for CSI AV. This is explained by the fact that 
the reference trajectory from which the IC's were derived rep- 
resented a 58 NM circular orbit, due to a smaller DOI burn than 
is currently being used. The difference between the DOI burn and 
the current value, is approximately the same as the difference 
between the CSI burn and its current value. 



3-U 


c < 

o 

»H F-1 
«P w 
crt CO 

& 

Q 

IS 

•3 H 

fi " 

-P 

CO 


OJ 

VO 


LA on v£) 

o av la 

• • • 

LA GO QO 


ao 


on 

on 


b- 

-3- 


* 

• • 

on vo 


oo 

LA 


CO 


CVJ 


-3 

t- 


OJ 


OJ 


vOOvvOrovoroHOv 

in aO OJ on ?▼> <J\ LA vO 


O 

•P 

0) 

W 


o 


w 

p 

a) 

I 

§ 


< 
Lj 1 

o 

<D 

CO 

& 

00 

LA 



B 

CO 

£ 

i 

CVJ 

rH 

-3* 

OJ 


OJ 

X 


CO 


OJ 

| 

X 



1 

on 

Cm 




O 



a) 




rH 

0) 



P 

'O 



as 

G 


o 

H 

■P 

< 

< v 


•H 

Eh 

CO 


I 

B 

CO 

£ 


X 

a> 



P 

<u 

S3 

e-s 

CO -p 


3 


c 

o 


G -P 
•H 2 

g H 

o o 

S CO 


-3- 

on 

3 


-3- 


-3- 

on 


t- 

-3- 


LA 

on 


o 

vO 


8 

on 

i 


LA 

CVJ 


& 

vo 


s 

LA 

CVJ 


VO 

-3 

OJ 


-* 

LA 


-4 

'C 


LA 


LA 


& 

on 


o 

OJ 

on 


LA 

LA 


vo 


R 


VO 


1/N 

rH 

• 

O 

LA 

• 

§ 

• 

8 

• 

rH 

ON 

• 

o 

CVJ 

• 

c- 

VO 

• 

on 

av 

• 

A- 

VO 

LA 

rH 

rH 

on 

_3 



OJ 


OJ 

os 

s 

So 

rH 

O 

00 

CVJ 

OS 

LA 

rH 

-3’ 


rH 


on 


CVJ 


$ 

CVJ 

I 


as 

trt 

TJ 

OJ 

a 

■p 

<p 

0 

co 

0) 

1 


OJ 

-p 

o 

co 

« 

<D 


O 

-P 

CO 

P 

<u 

<W 

<u 

p 

a> 


a 

$ 


<u 

-p 

O) 


*0 

0 

•H 

•P 

1 
-P 
•H 

<u 

.c 

•p 

s 

•H 

> 

•H 

T* 

>> 

P 







P <U -P 

■P GO) 

•H CO 
aS 

O -P QJ 

-3 

LA 

Os 

LA 

Os 

P P 

• 

t>- 

• 

rH 

• 

LA 

• 

on 

• 

§ 

CVJ 

-3 

-3 

C- 

OJ 

1 S3 

<D a co 

P > a> 


-P > 

O 4h 
-P o 

(0 P 
P <U 
<u P 
Ch a 
<d 3 
P G 


<U 
P 
<D -P 

* » 
p 


§ 

<u 

22 

fi 


0) CO 
■P 0) 

if 


CO 


CO 


co 




CO 


Q 


Q 


P 

ffi 



Q 


P 

S3 

P 

s 

<u 

> 

* 

P 


> 

> 

> 

> 

> 

> 

> 

> 

> 

> 

<i 

<J 

<3 

<3 

S 

<3 

<3 


<3 

<3 

rH 

rH 

OJ 

OJ 

G 

M 

5S 

K 

M 

M 

O 

O 

O 

o 

o 

o 

aJ 

CO 

Q 

P 

CP 

CP 

O 

a 

5S 

O 

U 

O 

Eh 

Eh 

X 

s 

3E 

as 



3-5 


c o 

o < l) 
•H CO 

P P -"s. 
cd o> P 
•H CO <M 
> 


o 

T3 Q 4> 
H to 
«9 P V 
'd <U P 
C CO <M 


8 

on 

oo 

,60 

S 

CVJ 

r— < 

s 

* 

s 



rH 


on 

ro 

CVJ 

CVJ 

H 

rH 


rH 

-4 

lA 

A 

rH 

C-* 

OV 

VO 

o 

CVJ 

on 

-4 

OV 

on 

IA 

rH 

o 

•H 


rH 







P 

(4 

- >» 

CO (0 


a h 

P o 

° t 

0> 0) 


P <u 

G 

0) o 

0 ) 

<y p 


o 

< 0 ) CVJ 

co o 

p \ 

fl) p 

CO <w l 

c 

cd 

0 ) 

S _ o 
Q 0) 
co 

P O 

<U P . 

CO <M O 


on vo 

on c— 


CT\ ON 

CVJ LTV 


C"- lA 

On CO 


00 o 
oo ia 


O co 

c b 

o 

> t 

.<1 Jm 

^ (II 


on 

<8 

O 

•H 

c 

1 

P 

o 

on 


o 

a> 

cd 

P 


rH 

O 

c 

P 


cd 

fl 

c 

a 

co 


a; 

P< 



O 


p 

Ph 


<8 

u 

w 


CO P 
<U 3 
O O 

c P 

<D P 

£ 
a) & 

<p 

<M 

•H 

Q 


ij 


o 

a) H- 

co 


O H lA VO 

-4 -4 VO rH 

CVJ CVJ CVJ CVJ 


lA- 4t—CJ\-4 

cvj on on irv on 


P* > 

> 

0 ) <3 

<3 

d ^ 

<3 

0 ) 


C M 

K 

« CO 

Q 

E 9 

O 


55 > 

> O 


CO 

2 : 3 25 

> > > 

<3 O <1 

< < < 


rH H CVJ CVJ 
M M U O O U 

fi h g g g U 


CO <u 

p 

g 1 

CO p 

0) o 

V. p 

& 

OJ 0) 

t« P 

P 

Q 

co 

P P 

<L> G 

CO <D 

CO 


tJ < 

a> 

p p 

C0 (U 
•H CO 
rH 

on a> 


cd p 

H co 


CO p 

H 3 


« 5 9 

P cd 

cd 'd *d 

^ • 8 

<d x: 

P P •* 


•H 
G »d 
a> a 
a> < v 

$ u 
V t3 

■° § 

a> 

o u 

c o 
<u co 
^ c 
a> <u 

dn CO 
€ 

•H P 


a) p 
p p 

+> s 

CO 

P 'd 

g <u 

<y -p 


0 ) 

Vi 

o 

> 

> 

<3 

<3 

<3 

M 

Hi 

CO 

O 

CO 

u 

O 

u 

fi 

o 

<u 

Cm 


oj 



> 



> 

oj 

p 

p 


3-6 


3.2 Miss Distance 

The miss distances were established by Sets B and C. 

The average in plane iuis3 distance at the point of closest 
approach for the 50 cases without errors (Set C) was 494 
feet, and for the 100 cases with errors (Set B) was 1664 feet. 

The average, mean, and standard deviation of the com- 
ponents of the miss distance in a local vertical coordinate 
system with X along the radius vector of the LM, Z along the 
angular momentum vector of the LM, and Y completing the right- 
handed system were as follows: 


Table 3-4 

Coordinates at Closest Approach 


Axis 


Average 


Mean 



SET B 

SET C 

SET B 

SET C 

SET : 


Feet 

Feet 

Feet 

Feet 

Feet 

X 

939 

194 

-329 

122 

1488 

Y 

1374 

454 

102 

392 

1952 

Z 

301 

251 

51 

17 

4ll 


Standard Deviation 
SET C 
Feet 


208 

256 


307 
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3.3 A V Used 

The mean total translation AV used in the 50 cases without 
errors (Set D) was 136.3 ft/sec with a standard deviation of 
15.0 ft/sec while the mean for the 100 cases with errors (Set A) 
was 150.5 ft/sec with a standard deviation of 21.0 ft/sec. 

The minimum and maximum AV cases without errors required 
120.5 ft/sec and 180.1 ft/sec, respectively, while with errors 
minimum and maximum AV cases required 120.2 ft/sec and 230.4 ft/sec. 
Figure 9 shows the distribution of total AV. 

A breakdown of how the AV was used is shown in the following 
table: 

Table 3-5 
Translation AV 


Maneuver AV (Set D) AV (Set A) 



Average Without Errors 
ft/sec 

Average With Errors 
ft/ sec 

CSI 

46.3 

47.1 

CDH 

13.6 

14.6 

TPI 

26.1 

27.5 

MCC1 

5.1 

8.9 

MCC2 

5.6 

9.9 

Braking and LOS 
Control 

39.7 

42.6 


Arrival Time at TPI 


The mean arrival time at TPI for the 50 cases without errors 
(Set D) was 23 seconds later than nominal with a standard 
deviation of 31 seconds, while the mean for the 100 cases with 
errors (Set A) was 21 seconds late with a standard deviation of 
4 minutes and 27 seconds. 

Figure 10 shows the distribution of arrival time at TPI 
over intervals of two minutes for the 100 cases with errors. 
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3.9 A T from TPI to TPF 

The mean AT from TPI to close approach without braking for 
the 50 cases in Set C was 2586 seconds with a standard deviation 
of 9 seconds while the mean AT of transfer without braking for 
the 100 cases in Set B was 2^84 seconds with a standard devia- 
tion of 44 seconds. 

The mean AT from TPI to intercept with braking and LOS 
control for the 50 cases in Set D was 2747 seconds with a standard 
deviation of 32 seconds while the mean AT of transfer with 
braking for the 100 cases in Set A was 2756 seconds with a standard 
deviation of 99 seconds. The nominal case required 2737 seconds 
with braking and LOS control. 

Figure 11 shows the distribution of the AT transfer with 
braking for the 100 cases with errors. 
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TERMINAL PHASE INITIATION 


FIGURE 3 
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F MISSION-RELATIVE REFERENCE TRAJECTORY 


FIGURE 6 
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